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Abstract 



We propose a new multiple integral representation for the correlation function 
(cfcTm+i) °f the XXZ spin-^ Heisenberg chain in the disordered regime. We 
show that for A = 1/2 the integrals can be separated and computed exactly. 
As an example we give the explicit results up to the lattice distance m — 8. It 
turns out that the answer is given as integer numbers divided by 2^ m+1 ^ . 
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In this paper we study the correlation functions of the XX Z spin-^ Heisenberg chain pQ. 
The Hamiltonian of this model is given by 



M 

H =J2 i a >™+l + «+! + A («+l - 1)} , (1) 
m=l 

where A is the anisotropy parameter and (Tm y ' z are local spin operators associated with the 
m-th site of the chain. The boundary conditions are periodical. 

In 2001 Razumov and Stroganov '% observed that, at the special value of the anisotropy 
parameter A = 1/2, the model exhibits a series of remarkable properties. It was conjectured in 
particular that, in the thermodynamic limit, the probability to find a ferromagnetic string of 
length m in the antiferromagnetic ground state is proportional to the number of alternating sign 
matrices of size mxm. This conjecture was proved in 3 using multiple integral representations 
of the correlation functions [I] , [3] , [H] , [7j . We use here the same approach to study the two- 
point correlation function of the third components of local spin. 

Following the strategy of [HI, 0, we introduce the operator 

Q^) = f[( 1 -^ + 1 -^-<), (2) 

n=l ^ ' 

where k is a complex number. The ground state expectation value of this operator is the 
generating function for the correlation functions of the third components of spin, namely, 



d 2 

= 2 -^D 2 m (Q K (m)) 



~ 1, (3) 

K=l 



where means the second lattice derivative. 

To study the generating function (Q K (m)) in the thermodynamic limit, we use the following 
multiple integral representation: 
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sinh(Afc — Zj — iQ 
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Here, 
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sinh(z - A) smh(z — A — iQ smh(z + i§) 
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and the integrals over the variables Zj are taken with respect to a closed contour T which 
surrounds the point — i(,/2 and does not contain any other singularities of the integrand. 

The equation (jlj) can be derived from the one obtained for {Q K {m)) in jjj by expanding the 
integrand into power series over k with forthcoming moving a part of the integration contours 
for variables Xj to the lower half-plane. 

The equation @ is valid for the homogeneous XX Z chain with arbitrary — 1 < A < 1. If 
we consider the inhomogeneous XX Z model with inhomogeneities £1, . . . ,£ m , then one should 
replace in the representation @ the function ip m in the following way: 



sinh(z - - i() 
sinh(z - Q 



. V (A) ^lisinh(A-^-<)- (6) 

6=1 6=1 v so s; 

In order to come back to the homogeneous case, one should set = —i(/2, k = 1, . . . , m in ©. 

In the inhomogeneous model, the integration contour T surrounds the points £i, . . . , £ m , and 
the integrals over ^ are therefore equal to the sum of the residues of the integrand in these 
simple poles. Computing these integrals and setting £ = 7r/3, we obtain an alternating sum 
over the permutations of the set {£i, . . . ,£ m }: 



(-1) 2 
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a>b 



sinh3(£ a -£ b ) 



sinb/ ! (£ a -£& 



?7T N 

3 / 



Y ^ m ~ n c^Y{-i) [p] TT 11 



n=0 



n smh(A a - a 6 ) n Kt, 

a>b j=l 



p(j)' ^7 



' n <(Ai,4p(,-)) 

j=n+l 



(7) 



]1 < 11 sinh(A a - £ p(6) 

a=l 6=1 



»7T ' 
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11 sinh(A a - £ p(6) + f ) 



6=n+l 



It is easy to see that the double products containing the parameters {A} can be written as a 
Cauchy determinant with sinh _1 (Aj — ^ v {k) — §■) in the first n columns and sinh _1 (Aj — i p (k) + ^f") 
in the remaining ones . Thus, the integrand can be presented as a determinant of a block-matrix: 



(<9 K (m)) 



sinh3(£a - 6) 




-"^En n 

n=0 p a=l 6=n+l 



sinh(^ p(f)) - £ p(a) - f ) sinh(£ p(a) - f p(6) ) 
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Here the sizes of the blocks in the determinant are respectively: nx n; nx (m — n); (to — n)xn; 
and (to — n) x (to — n). 

We see that the original m-fold integral over Xj is now factorized. Indeed, Xj enters only 
the j-th line of the determinant and, hence, one can integrate each line separately. Thus, we 
arrive at 

«■("» - w n TJlt "&) t " m - e det4<»)(K 7+ }, Kl _}) 

a>6 Smn ^ a W n=0 «}={«7+}u{«7_} 

7 + |=n 



» n n 



S mh(6-^-f)sinh(^ a -6) 



sinh 2 (£ b - £ a + 
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with 
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Here the sum is taken with respect to all partitions of the set {£} into two disjoint subsets 
{£7+ } U {£ 7 _ } of cardinality n and m — n respectively. The first n lines and columns of the 
matrix <I>( n ) are associated with the parameters £ € {£7+}- The remaining lines and columns 
are associated with £ 6 {£7-}- 

Thus, we have obtained an explicit answer for the generating function (Q K (m)) of the in- 
homogeneous XXZ model. However, the homogeneous limit of @ is not obvious, since each 
term of the sum in Q becomes singular in the limit 1 £*, — * 0, k = 1, . . . , to. In order to check 
that the sum over partitions in the r.h.s. of Q remains indeed finite in the limit £fc — ► 0, one 
can introduce again a set of auxiliary contour integrals: 
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sinh(z,; 
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n-^r sinh(z a - z b - f ) sinh(z a - z b + . 
II sin JL ^ 5_-det$(^-^). (11) 

0=1 b=n+l Smh ^ a - Z b) 



Since (JSJ depends only on the difference — one can obtain the homogeneous XXZ chain in the limit 
£fc — > as well as in the limit ^ — > — i£/2. 
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Here the integration contours T{£ =p ^} surround the points {£ — ^} for zi, . . . , z n and {£ + ^} 
for z n +i, ... ,z m respectively. It is easy to check that, taking the residues of the integrand in 
the simple poles {£=t ^}, we immediately reproduce ©• On the other hand, one should simply 
set £jfc = 0, k = 1, . . . , m to proceed to the homogeneous limit in 1|11|) . As a result we obtain 
poles of order m in the r.h.s. of 1)11(1 . 

Certainly, the remaining integral is of Cauchy type and, after the change of variables xj = 
e 2z j ( it reduces to the derivatives of order m — 1 with respect to each Xj at x\ = ■ ■ ■ = x n = e~ 
and x n+ \ = ■ ■ ■ = x m = e~~ . If the lattice distance m is not too large, the representations @, 
(fTTj) can be successfully used to compute (Q K (m)} explicitely. As an example we give below the 

2 

list of results for P„i(k) = 2 m (Q K (m)) up to m = 9: 
Pl(«) = l + K, 

P 2 ( K ) = 2 + 12k + 2k 2 , 

P 3 ( K ) = 7 + 249k + 249k 2 + 7k 3 , 

P 4 {k) = 42 + 10004k + 45444k 2 + 10004k 3 + 42k 4 , 

P b ( k ) = 429 + 738174k + 16038613k 2 + 16038613k 3 + 738174k 4 + 429k 5 , 
P 6 ( K ) = 7436 + 96289380k + 11424474588k 2 + 45677933928k 3 + 11424474588k 4 
+96289380k 5 + 7436k 6 , 

P 7 {k) = 218348 + 21798199390k + 15663567546585k 2 + 265789610746333k 3 

(12) 

+265789610746333K 4 + 15663567546585k 5 + 21798199390k 6 + 218348k 7 , 
P 8 (k) = 10850216 + 8485108350684k + 39461894378292782k 2 
+3224112384882251896k 3 + 11919578544950060460k 4 + 3224112384882251896k 5 
+39461894378292782k 6 + 8485108350684k 7 + 10850216k 8 
P 9 (k) = 911835460 + 5649499685353257k + 177662495637443158524k 2 
+77990624578576910368767k 3 + 1130757526890914223990168k 4 
+1130757526890914223990168k 5 + 77990624578576910368767k 6 
+177662495637443158524k 7 + 5649499685353257k 8 + 911835460k 9 . 
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Using the results (fl2|) . we can calculate the correlation function up to m = 8: 





= -2~\ 




= 7 - 2" 6 , 


Kd> 


= -401 • 2~ 12 , 




= 184453 • 2~ 22 , 




= -95214949 • 2~ 31 , 




= 1758750082939 • 2" 46 , 




= -30283610739677093 • 2~ 60 , 




= 5020218849740515343761 • 2 



(13) 



For m = 2, 3 we reproduce the results (Hj, JOJ where the two-point correlation functions of the 
XXZ chain were calculated up to the lattice distance m = 3 for general A. 

One can also compare the exact results (|13j) with the values given by the asymptotic formula 
suggested in [TT] . 



m 


Exact 


( a i a m+i) Asymptotics 


1 


-0.5000000000 


-0.5805187860 


2 


0.1093750000 


0.1135152692 


3 


-0.0979003906 


-0.0993588501 


4 


0.0439770222 


0.0440682654 


5 


-0.0443379157 


-0.0444087865 


6 


0.0249933420 


0.0249365346 


7 


-0.0262668452 


-0.0262404925 


8 


0.0166105110 


0.0165641239 



One can see that the correlation function (ofc>"m+i) approaches its asymptotic regime very 
quickly and that, already for m = 4, the relative precision is better than 1%. 

In conclusion, we would like to comment on the results obtained in this paper. We have 
shown that, for the special value of the anisotropy parameter A = 1/2, the multiple integral 
giving the generating function {Q K {m)) can be factorized and computed explicitly in terms of 
the inhomogeneities. It was shown recently in the papers ^2], ^H] that a factorization of the 
multiple integral representations for the correlation functions of the inhomogeneous XXZ chain 
is possible for generic value of A. However, as in ^3] , we see that it is not straightforward to 
take the explicit homogeneous limit of @ for general m, although it is possible to do it for 
small distances. One way to obtain explicitly this homogeneous limit is to come back to the 
multiple integral of Cauchy type (|11[). but with non- factorized integrand. 
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Concerning the exact numerical formulas we obtained for the generating function (Q K (m)) 
up to ra = 9, we would like to mention the remarkable fact that all the coefficients of the 
polynomials P m («) turned out to be integers. We conjecture that this property holds for all 
m. In particular, the highest and lowest coefficients of P m {^) correspond to the emptiness 
formation probability and, hence, are equal to the number of alternating sign matrices of the 
size mx m 0, [5J. For the other coefficients, this property should be related to the structure of 
the ground state at A = 1/2 conjectured by Razumov and Stroganov in [2] (see also [Tl], [To] . 
|16|). It would be very interesting to reveal the combinatorial nature of these coefficients. 
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